##### SCRIPT QUE GERA A TABELA 6 #####

# 1. CARREGA PACOTES ------------------------------------------------------

library(tidyverse)
library(janitor)
library(naniar)
library(mice)
library(xlsx)

# 2. CARREGA DADOS ------------------------------------------------------

#' Nesta etapa, foi implementada a leitura dos dados combinados de Enade-RAIS

# 3. IMPLEMENTA MODELOS DE REGRESSAO LOGISTICA --------------------------

# 3.1 Modelo geral (sem separar por areas de estudo) ----

model_emp_0 <- with(dados_graduados_empregados, glm(ind_empregado ~ 
                                                          1,
                                                    family = binomial))

model_emp_1 <- with(dados_graduados_empregados, glm(ind_empregado ~ 
                                                          tp_sexo + 
                                                          raca_cor + 
                                                          idade_faixa + 
                                                          edu_fam,
                                                    family = binomial))

model_emp_2 <- with(dados_graduados_empregados, glm(ind_empregado ~ 
                                                          tp_sexo + 
                                                          raca_cor + 
                                                          idade_faixa + 
                                                          edu_fam + 
                                                          nt_ce_quartil + 
                                                          nt_fg_quartil +
                                                          trab_grad,
                                                    family = binomial))

model_emp_3 <- with(dados_graduados_empregados, glm(ind_empregado ~ 
                                                          tp_sexo + 
                                                          raca_cor + 
                                                          idade_faixa + 
                                                          edu_fam + 
                                                          nt_ce_quartil_area + 
                                                          nt_fg_quartil +
                                                          trab_grad + 
                                                          igc_faixa_4c +
                                                          setor_ies_bin + 
                                                          regiao_curso,
                                                    family = binomial))

# 4. GERA TABELAS DE RESULTADOS -------------------------------------------

combine_emp_0 <- pool(model_emp_0)
summary(combine_emp_0)
combine_emp_0_fit <- combine_emp_0$glanced

combine_emp_1 <- pool(model_emp_1)
summary(combine_emp_1)
combine_emp_1_fit <- combine_emp_1$glanced

combine_emp_2 <- pool(model_emp_2)
summary(combine_emp_2)
combine_emp_2_fit <- combine_emp_2$glanced

combine_emp_3 <- pool(model_emp_3)
summary(combine_emp_3)
combine_emp_3_fit <- combine_emp_3$glanced

combine_emp_3_out <- data.frame(summary(combine_emp_3) %>% 
                                   mutate(estimate = exp(estimate)))

# 5. EXPORTA TABELA -------------------------------------------------------

xlsx::write.xlsx(combine_emp_3_out, 
                 file = "./output/Tabela-6.xlsx", 
                 sheetName = "raw", 
                 row.names = F, 
                 append = T)
